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INTRODUCTION 


AD-P003  648 


The  subject  matter  of  the  first  part  or  mis  study  is  the  deri- 
vation of  coupled  bending-torsion  relations  characterizing  the  dyna- 
mical behaviour  of  unsymetrical  cross-section  beams.  This  allows  for 
the  further  definition  of  a beam  element  with  seven  degrees  of  free- 
dom per  node.  The  numerical  results  obtained  with  the  FEM  are  compa- 
red to  experiment  in  some  section  shape  cases. 

In  order  to  characterize  the  displacement  field  of  considered 
beams,  the  method  of  integrated  displacements  allows  us  to  consider 
the  so-called  secondary  effects!  longitudinal  warping  inertia,  and 
shear  deformation  due  to  both  shearing  forces  and  nonuniform  warping. 
The  literature  on  dynamical  flexure  and  torsion  of  beams  is  extensi- 
ve. Cowper  •'£-4-^  for  flexure  and  Gay-Boudet  *£2^for  torsion  introduce 
integrated  displacements.  We  extendin  the  following  this  notion  to 
the  study  of  coupled  bending- tors i on  of  an  homogeneous  straight  beam. 
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BASIC  THEORY 


f 


> t 


l • 


i f 


2 . 1 Displacement  Field 
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(2  is  the  domain  occupied  by  the  cross-section.  T is  the  bounda- 
ry, and  G is  the  centroid.  In  the  plane  of  the  section,  the  principal 
axis  are  noted  Gx2»Gx3.  C (02,03)  is  the  shear  center,  as  defined  by 
Trefftz  [3] (figure  1).  ^ 

The  displacement  of  any  point  M of  the  beam  is  noted  XM(x^) , £ • 

with  components  X£  (i-1,2,3).  Let  us  define  the  seven  displacement 
parameters : 


Three  angular  parameters: 

8i(xi  , t ) - 7 / (GMxXm)  . X£  dS2 

li  a 
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Three  linear  parameters: 

Ui(xi , t)  = ^ /xM.xi  dQ 
b a 

d-b) 

A warping  parameter: 

§ (x 1 , t)  = — /$ . x 1 da 

I$  a 

d-c) 

t 
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$(x2,X3)  is  the  Saint-Venant  warping  function  defined  in  C,  and 
I$  = JV  da  is  the  quadratic  warping  moment. 
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x = r + gmx^  + $§xx  + n 
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According  to  the  above  definitions,  necessary  orthogonality  conditi- 
ons hold  in  the  domain  ft  for  the  functions  1,  X2 , X3,  $,  n 1 » together 
with : 

/ n2  dft  » /n 3 dft  - /(x2 n3  - x3n2)  dft  - o (3) 

q ft  ft 


2 . 2 Equations  of  Motion 

Assuming  the  lateral  surface  of  the  beam  free  of  any  force,  we 
take  into  account  as  first  hypothesis  (H . 1 ) the  assumption  that  nor- 
mal stresses  022  and  033  are  expected  to  be  negligible  compared  to 
Oil.  Then,  the  principle  of  virtual  work  associated  with  the  displa- 
cement field  (2)  leads  to  the  classical  set  of  motion  equations: 
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F.  1 + P 


(4-a,b, c) 


[l]$> tc  " | + m + XjXF 


(5-a,b , c) 


And  the  seventh  motion  equation  is  the  bimoment  one; 

P I , t t **  ®*1  + k + (Mj—  GJ0j,j)  — *-2^3  ^ ^3^2  (^) 

in  which  B ■ / $o  l j dft  is  the  generalized  bimoment, 
ft 

2 . 3 Constitutive  Equations 

In  order  to  obtain  a technical  formulation  for  the  constitutive 
equations,  we  assume  as  second  hypothesis  (H,2)  that  the  change  in 
the  deformation  ri  of  two  infinitely  near  adjacent  cross  sections  is 
neglected  in  order  to  evaluate  the  longitudinal  and  shear  stresses 
a 1 j , j » 1,2,3.  Thus , starting  from  the  Hooke's  law  for  a linear  elas- 
tic body,  the  constitutive  equations  below  are  deduced  from  integra- 
tions over  the  domain  of  the  section: 
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(7-a,b,c,d) 


Moreover,  noting  that  the  first  component  Hi  of  r|  (2)  is  essentially 
due  to  the  shear  forces  F2 , F3  and  nonuniform  warping  moment  noted 
Mnuw  * (Mi-GJ9i,i),  we  obtain  the  three  followed  coupled  relations 
written  in  matrix  form: 
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KX1  -Kl2  -K13  1 r Mnuw/  G ( I 1 — J ) 


Ui,  1 + 0,+  c,0 
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F2/GS 


F3  / GS 


(8-a,b , c) 
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In  (8)  appear  nine  shear  coefficients  with  the  symmetry  proper- 
ties : 
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Kij  = Kjl  » 3=2,3  (9-a,b,c) 


Starting  from  the  local  equilibrium  of  continuous  media,  we 
show  that  the  shear  coefficients  K£j  can  be  computed  after  solving 
three  Poisson’s  problems  in  the  cross  section,  namely: 


V2g  ■ f(x2,x3)  over  fi,  with  successively  f « x2,x3,$, 

3g/9n  ■ 0 along  the  boundary  T (10-a,b,c) 

with  j g dft  ■ 0 
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! For  the  complete  expressions  of  K^j,  see  appendix  A.  The  effective 
' calculus  of  all  bending-tors  ion  constants  involved  in  the  aforemen- 
tioned relations  ( X £ , J,  1$,  K^j , C£>  has  been  performed  by  means  of 
a boundary  element  method  [A]  . 

3.  FINITE  ELEMENT  FORMULATION 

For  the  setting  up  of  a finite  element  formulation,  let  us 
first  look  at  the  technical  expressions  of  potential  and  kinetic 
energies . 

3 . 1 Poten t i a 1 S tr ein  Energy  V(xi , t) 

For  a beam  element  of  length  dxi,  and  according  to  H.l,  the 
| general  form  of  strain  energy: 

dV  - { j Joi jEijdnJdxj 

n 

reduces  to: 

; V,,  - - /(One  u + 2012e12  + 2ol3el3)  OD 
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Considerable:'  of  (2)  and  H.2  leads  to  the  practical  expressions: 
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v.i  - V0M  + vT, 
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in  which: 
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( 1 3-b) 


3 . 2 Kinetic  Energy 

In  the  same  way,  we  retain  the  simple  expression: 

T,,  - \ pS(tj,t)2+  \ pl*$?t  + \ p{^,t}T|l|  {S,tJ  (14) 

3 . 3 Definition  of  the  Element 

Seven  displacements  per  node  k characterize  the  motion  of  ti;. 
two-nodes  beam  element.  In  a matricial  form:  , 

{Wk}  - [u1,u2,e3,u3,02,01,5]T 


Substituting  the  static  case  for  the  quasi-static  dynamic  one 
leads  to  an  interpolation  matrix  [a]  not  detailed  here,  allowir  a 
displacement  field  in  the  form  of: 


{ W (x j , t) } - [ A ( x J ) ] { W ( t ) > 


(15) 
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3 . 4 Stiffness  Matrix 

Introducing  the  displacement  field  (15),  the  element  strain 
energy  is  obtained  after  integration  of  (13)  over  the  length  l of  the 
element,  and  written  in  the  matrix  form: 

2V  - {W}T[d]{W} 

We  detail  in  Appendix  B the  stiffness  matrix  [d]  in  the  simpli- 
fied case  i^j+Kjj'vO.  ’it  can  be  noted  that,  for  symetrical  section  sha- 
pe cases  and  for  uniform  warping  (Saint-Venant  torsion),  this  stif- 
fness matrix  reduces  to  the  classical  one  derived  by  Przemieniecki  [5J . 

3 .5  Mass  Matrix 

The  calculus  of  the  kinetic  energy  for  the  whole  element  by 
means  of  (15)  leads  to  the  form: 


2T  = {W , t }T[m]  (W  , t } 


in  which  [m]  is  a consistent  mass  matrix,  reached  after  very  heavy 
calculus.  In  the  much  more  simpler  case  of  lumped  mass  approximation, 
the  kinetic  energy  reduces  to: 
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from  which  a simple  diagonal  mass-matrix  is  easily  derived. 

4.  APPLICATIONS 

Numerical  and  experimental  tests  have  been  performed  on  several 
cantilever  beams.  For  the  dynamical  flexure  and  torsion  of  a rectan- 
gular section  beam,  numerical  values  of  frequencies  provided  by  FEM 
and  analytical  ones  are  in  good  agreement  with  the  experiment,  espe- 
cially for  high  frequencies.  The  dynamical  torsion  of  aal-section 
beam  has  been  likewise  investigated.  In  such  a case,  exact  analytical 
values  of  torsional  constants  cannot  be  reached,  needing  a previous 
computational  work.  The  results  for  dynamical  case  are  detailed  in 
reference  [6],  and  some  of  them  are  recalled  in  appendix  C.  In  a same 
manner,  we  have  also  tested  two  U-section  beams  (thick  and  thin).  Here 
coupling  between  flexure  and  torsion  occurs,  and  the  whole  theory  abo- 
ve applies-  We  shall  present  both  numerical  and  experimental  frequen- 
cies. Rather  good  agreement  can  oe  noted  for  all  tests, the  rank  of 
the  modes  concerned  depending  on  the  number  of  elements  of  the  discre- 
tization. 

5.  CONCLUSION 

The  study  of  coupled  bending-torsion  can  be  performed  with  ac- 
ceptable accuracy  by  means  of  the  formulation  above,  starting  from 
the  definition  of  seven  displacements  parameters  in  each  section  of 
the  beam.  The  finite  element  derived  allows  a simple  numerical  pre- 
diction of  the  dynamical  behaviour  of  beams  with  any  cross  sections. 
Nevertheless,  wc  must  keep  in  mind  that  an  accurate  computation  of 
coupled  bending-tor sion  constants  is  the  first  stage  when  using  this 
e lemant . 


NOMENCLATURE 

E,G  Young's  and  shear  modulus. 

S,I,,I2,I3  Cross  section  area  and  quadratic  moments  of  inertia  in  G. 
I ,j,  Quadratic  warping  moment. 

J Saint-Venant ' 8 torsional  rigidity. 

Normal  and  shear  forces. 

MlfM2,M3  Torsion  torque  and  bending  moments. 

U,$  Linear  and  angular  displacement  vectors  (components : u^  and  0^) . 
x^  Unit  vectors  of  principal  axis, 

p Mass  per  unit  volume. 

( ),£  Partial  derivative  with  respect  to  x£. 

( ) , t id . with  time  t . 

3 / 3n  Outward  normal  derivative  along  the  boundary  F. 

p,m,b  Distributed  forces,  moments  and  bimoment  along  the  beam. 

[i]  Diagonal  inertia  matrix  of  cross  section, 

rf  Complementary  displacement  vectc.% 
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/ppendix  A 

Shear  Coefficients  K^j 

Solutions  of  (10-a,b,c)  allow  the  knowledge  of  the  functions 
gi0(i*l,2,3)  verifying  over  the  domain  0 of  the  cross-section: 
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<J>  - — x. 


v2g3 
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with  3gi0/3n  = 0 ; i = 1 ,2,3  along  the  boundary  T. 
Then,  the  K^j  are  deduced  in  the  form  of: 

K 1 1 = ~~  /$ g 3 0 ^ 


and  for  i,j=2,3  ; i^j  ; k=j-l  : 
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The  properties  of  symetry  of  the  k£j  are  shown  by  means  of  the 
Green-Gauss  theorem,  and  we  find: 
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Appendix  B 

Stiffness  Matrix  [d] 

He  introduce  the  following  dimensionless  notations: 


,2,2  GJ4‘ 

0 ' 1 * K"  1 x t ' oil*  » k.  * 
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Then,  the  symetry  properties  of  [d]  leads  to  the  following 
terms,  in  the  simplified  case  : ij*  j^Kj  j "0 . 
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Appendix  C 


The  study  of  natural  torsional  frequencies  of  a cantilever  beam 
performed  by  means  of  the  finite  element  procedure  has  been  compared 
to  the  experimental  datas.  The  results  concerning  the  relative  error 
are  shown  on  table  1 below. 
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! MODE  N* 

1 

2 

3 

4 

5 

i 

7 

8 

■ 

D 

ERROR  X 

0. 

0.9 

0 . 

0.7 

-1 . 

-0.9 

-1  .6 

-1.3 

0.1 

2 . 

Table  1. 
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